Long life of Gauss-Bonnet corrected black holes 
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Dictated by the string theory and various higher dimensional scenarios, black holes in D > 4- 
dimensional space-times must have higher curvature corrections. The first and dominant term is 
quadratic in curvature, and called the Gauss-Bonnet (GB) term. We shall show that although 
the Gauss-Bonnet correction changes black hole's geometry only softly, the emission of gravitons is 
suppressed by many orders even at quite small values of the GB coupling. The huge suppression of 
the graviton emission is due to the multiplication of the two effects: the quick cooling of the black 
hole when one turns on the GB coupling and the exponential decreasing of the grey-body factor 
of the tensor type of gravitons at small and moderate energies. At higher D the tensor gravitons 
emission is dominant, so that the overall lifetime of black holes with Gauss-Bonnet corrections is 
many orders larger than was expected. This effect should be relevant for the future experiments at 
the Large Hadron Collider (LHC). 
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I. INTRODUCTION 

During the past decade high energy physics received 
a great impact from theories implying existence of extra 
dimensions in the world. These are the string theory [1] 
and higher dimensional brane- world scenarios [2]. The 
low energy limit of the string theory can be described 
by the slope expansion in powers of the inverse string 
tension (or of the inverse square of the fundamental string 
scale that produces higher curvature corrections to 
the Einstein action. The quadratic term in curvature 
(given by the so-called Gauss-Bonnet invariant) is the 
leading correction that can affect the graviton excitation 
spectrum near the flat space. 

The extra dimensional scenarios also suggest that the 
fundamental gravity scale might be around the weak 
scale ~ TeV . Thus, at particle collisions with the cross 
section ^ Trr^, where is the Schwarzschild radius, and 
energies larger than Af,, the production of mini-black 
holes should start. These black holes are intrinsically 
higher dimensional and usually modeled by the Tangher- 
lini metric, which is the solution of the D-dimensional 
Einstein equations. However, in order to have a mathe- 
matically noncontradictory gravity in higher dimensions, 
one has to take account of higher curvature corrections of 
the same form as those appearing in the slope expansion 
of the string theory. The spherically symmetric solution 
describing neutral static black holes in the D-dimensional 
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Einstein gravity with the GB corrections was obtained 
in [3]. This solution contains small corrections to the 
D-dimensional Schwarzschild- Tangherlini geometry and 
consequently properties of such Gauss-Bonnet corrected 
black holes were expected to differ only slightly from the 
Schwarzschild's ones. This happens, for instance, for the 
spectrum of proper oscillations of these black holes [4]. 

Unlike astrophysical black holes, whose Hawking evap- 
oration is negligibly small, mini-black holes are inten- 
sively evaporating what leads to the very short lifetime 
of these black holes, once they are created. The lat- 
ter is due-to strong production of various particles from 
the vacuum around a black hole and emission of them 
through the mechanism of Hawking radiation [5]. At 
large number of space-time dimensions D, the specific 
"tensorial" type of gravitons (respectively the D — 2 ro- 
tation group) dominates in the emission process [13]. 
Up to now, an impressively extensive literature is de- 
voted to the calculations of Hawking evaporation of the 
Schwarzschild- Tangherlini and Myers-Perry black holes 
[6-10], while evaporation of their higher curvature cor- 
rected generalizations was touched upon only in a couple 
of works [11, 12]. In particular, T. Rizzo estimated the 
energy-emission rate for the higher curvature corrected 
black holes, assuming that the grey-body factor equals 
unity [11]. This was expected to give the correct answer 
about the order of the intensity of the Hawking emission. 
Though, as we shall show in this paper, the contribution 
due to the grey-body factors can also considerably change 
the results. In [12], the scattering of Standard Model par- 
ticles around Gauss-Bonnet black holes was considered, 
though the calculations were terminated at the grey-body 
factors and the numbers of particles per frequency. Thus, 
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none of the above works calculated the energy-emission 
rate for Gauss-Bonnet black holes that is necessary for 
the estimation of the total emission of energy and thus of 
the black hole lifetime. Here we shall fill this gap and cal- 
culate the energy-emission rates for fields of various spin, 
including gravitons, and thus shall estimate the lifetime 
of Gauss-Bonnet black holes. 

In this work, we shall show that due to a number of 
reasons, the emission of the tensor type of gravitons is 
greatly (in fact exponentially) suppressed when one turns 
on the GB coupling a'. Thus, even at small values of the 
GB coupling constant a' the graviton emission is sup- 
pressed by many orders. This means that small GB cor- 
rections lead to a much longer life of higher dimensional 
black holes than was expected [13-15]. At first sight this 
enormous suppression would not seem trustworthy: why 
do slight corrections of geometry produce a very strong 
effect on the evaporation process? The reasons for this 
are "multiplication" of the two factors. First, the black 
hole gets much colder when one turns on the GB coupling 
and the emission rate is quadratic in temperature. Sec- 



ds^ = f(r)de - 



ond, the emission is proportional to the grey-body factor 
which is exponentially suppressed for tcnsorial gravitons. 
This explanation certainly did not make us trust the re- 
sult immediately. Therefore we reproduced our accurate 
numerical calculations by the semianalytical WKB esti- 
mations. 

The paper is organized as follows: Sec II briefly dis- 
cuss the deduction of wave equations for perturbations 
of fields of various spin. Sec III is devoted to numerical 
calculation of the coefficients of transmission, while Sec. 
IIIB gives WKB values of the coefficients. In Sec. IV the 
obtained scattering data are used for the calculations of 
the energy-emission rates. Using the WKB arguments 
Sec. V explains why the found enormous suppression of 
Hawking evaporation occurs. In Sec. VI we estimate 
the lifetime of Gauss-Bonnet corrected black holes and 
outline the future perspective for this direction. 

We shall consider the canonical ensemble, which leads 

to the same results as the microcanonical one if the black 
hole mass M is at least a few times larger than M* [11]. 



(1) 



- r''dnl_^, (2) 



I 

II. THE WAVE LIKE EQUATIONS 

A. Gravitational perturbations 

The Einstein-Gauss-Bonnet action is 

/ = -^^^ j d^'x^gR + a' J d'' x^{RabcdR'"""' - ^RcdR'" + R^), 
where a' is the coupling constant; a' = 1/2-k(.1. The metric has the form, 



^^"^'- ^ a{D-2,){D-Ay «^^)-yi+ (zp_2)r^-i ^"^^ 

where dVL\,_2 is the line element of a unit [D — 2)-sphere and a = IQirGoct' ■ As we mentioned earlier, we shall 
measure all quantities in terms of the black hole horizon radius, which we shall denote as ro, so that the black hole 
mass can be written as 
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The linearized perturbations of (2) can be written in the general form as 

gab gab + 5 gab = Qab + hab- (5) 

Indices of hab are raised using the background metric; therefore 
The first order variation of the Riemann tensors is: 

SRab"'^ = \ [Rab^^hf" - Rab^h/' + (VfaV^/l/ - VaV^/lfe'^) + a'^'^K' - ^b^"ha') } . (6) 
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FIG. 1: Square root of the reflection coefiicient \Zo/Zi\ for D = 6, a = 2, 1 = 2 tensor-type gravitational perturbations. The 
left panel shows the coefficient calculated by fitting the numerically solved equation (blue) and using the 6-th order WKB 
formula (red). In the right panel we plot the difference between the coefficients calculated using these two methods. 




FIG. 2: The effective potential V{r) (left panel) and the square root of the reflection coefficient \Zo/Zi\ (right panel) for D = 6, 
a = 2, I = 2 scalar-type gravitational perturbations. 



Using the symmetry of the background metric, one can decompose the perturbation equations into scalars, vectors 
and tensors according to the rotation group on the {D — 2)-sphere [20]. Then, the separation of angular variables is 
possible for this case [20]. 



After separation of the angular variables and implying stationarity, ^(t,r) 
gravitational perturbations can be reduced to the wavelike equations [20] , 



e *"*$(r), the dynamics of the 



— +c^^-y,(r) ) $,(r) =0, dn = 



dr 

W) 



(7) 



with the effective potentials which have a very cumbersome form. After some algebra, we managed to simplify the 
potentials obtained by Dotti and collaborators in [20] for the tensor, vector, and scalar types of the gravitational 
perturbation, respectively: 
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VT{r) = I{r)\ U - + , , , ^ JrD-^A{r)q{r), (8) 

r \ A{r) J sJr'^-^A{r)q(r) dri V 

Vv{r) = f{r)^-£^A{r) + ^r^-^A{r)qir)-^ / 2^.. ^ . ^ ' (9) 

V ,yr^-^A{r)q{r) 

^ 64r2(£)-3)2A(r)2g(r)8(4cg(r) + (£)-l)ii(g(r)2-l))2 ^ > 

Here we used the following dimensionless quantities: 

1/1 1 \ /I 1 



A{r) = 

B{r) = A(r)^(l + ^^) + {l 
R 



2 A I \ 1 



D-Aj \ D-4 

2 



.2^ 



a(L>-3)(£)-4)' 

Z7(r) = 5(L'-l)*5i?2(i_^^^_3(£,_1^5^((£^_-^^^2_^24c(l + i?))g(r) + 

+2(D - l)^(24c(£) - l)i?2 + 168c2(l + ii) - (£» - l)i?^(-3 + 5i? + 7£»(1 + R)))q{r) 

+2(D - l)''i?(-184c2 + iD- 1)(13 + L')i?2 + c(-84 + 44i? + 84D(1 + R)))qirf + 
+{D - l)3(384c3 - 48c(2 + D{3D - b))R'^ + 192c2(-ll + D+ (-15 + D)R) + 
+{D - l)i?2(-3(7 + 55i?) + £1(26 + 106i? + 7D{1 + R))))q{r)^ + 
+{D - lfR{-Mc^{D - 38) + {D- 1)(71 + D{7D - <dQ))R^ + 
+16c(303 + 255i? + 12,D'^{1 + R) - 2D{73 + 81R)))qirf + 
+4{D - l)2(96c3(-7 + D)- 8c{D - 1)(145 - 74D + 6D'^)R'^ - 
-8c^{9 - 175R + D{-58 - MR + IID{1 + R))) + {D - l)R'^{-5{79 + 2SR) + 
+L>(5(57 + 41i?) + i:>(-81 - 89R + 7D{1 + R)))))q{rf - 

-4(£) - l)2i?(8c2(43 + (72 - 13Li)D) + {D - l)(-63 + D(99 + D{-49 + 5D)))R^ + 

+4c(321 + 465i? + £>(121 - 39i? + £)(-123 - 107i? + 17D{1 + E)))))g(r)^ + 

+{D - l)(128c3(-9 + D){D - 5) + 32c(D - 1)(246 + D{9 + D(-55 + 8D)))R'^ + 

+6Ac^{D -5)(D'^-3 + (49 + {D - 4)D)R) - {D - l)R^{n73 + 565R + 

+L»(-4(997 + 349i?) + £1(6(393 + 217i?) + L»(-548 - 452i? + 45L>(1 + R))))))q{r)^ + 

+{D - l)R{-64c^{D - 5)(36 + D(-13 + 3D)) + 

+{D - 1)(635 + £)(-1204 + 3£)(294 + D{-92 + 9D))))R^ - 

-8c{D - 5)(63 + 31i? + D{127 + 191R + £)(-47 + D + (-79 + D)R))))q{rf + 

+2{D - 5){64c^{D -5){D-3) + 8c{D - l)(-27 + I?(141 + (-43 + D)D))R^ + 

+8c^{D - 5)(-3 + 77i? + D{D -2 + {D- 18)R)) + {D - 1)^R^{-SS{R - 7) + 

+L»(59 + 43i? + £'(-59 - 35i? + 9L»(1 + R)))))q{r)^° - 

-2{D - 5)i?(24c2(-ll + D){D - 5){D - 3) + (£> - l)2(-65 + £>(81 + D{7D - 39)))i?2 + 

+12c(-7 + D){D- 5){D - 3)(D - 1)(1 + i?))g(r)" + 
+ {D - 5)2(-l + D)R'^{mc{26 +{D- 9)D) + 
+ {D - 1)(77 -3R + £(-18 + D + {D- 2)R)))q{rf^ + 
+{D - 5)^{D - 3)2(£ - l fR^q{rY^, 

X = {D — 2){c + 1) = i{i + D — 3) is the eigenvalue of the angular part of the Laplacian, £ = 2, 3, 4 

A wide class of static black holes in the Einstein gravity is known to be stable against gravitational perturbations 

for any number of space-time dimensions [19]. This is not true for black holes in the Einstcin-Gauss-Bonnct theory 
which suffer from instability for sufficiently large (in units of the radius of the event horizon) values of the GB coupling 
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a [20, 21]. Wc shall consider here only small enough values of the GB coupling a which are below the threshold of 
the gravitational instability [20, 22]. In other words, we shall be limited by not very small masses of black holes M, 
which usually are at least one order larger than M,. 



B. Test scalar field 

We shall also consider the test scalar field in the background (2) which satisfies the Klein-Gordon equation 

V"Va$ = 0. (11) 

After separation of the angular variables the equation (11) can be reduced to the wavelike form (7) with the following 
effective potential 

VM^^ = + '-^^^fir} + ^/'m) . (-) 

where X = i{i + D ~ 3), i ^ 0, 1, 2,3,.. .. 

Let us note, that unlike the Schwarzschild black hole, the Gauss-Bonnet black hole has the effective potential for 
the scalar field (12) which differs from the effective potentials for the tensor-type gravitational perturbations (8). In 
the limit a — >■ (12) coincides with (8). 



C. Brane-localized fields 



In addition, we shall consider the Standard Model fields (scalars, fermions and gauge bosons) living on a 4- 
dimensional brane, which is embedded in the background of the Gauss-Bonnet black hole. The induced metric 
on the brane is given by a projection of the metric (2) onto the 4-brane [23] 

ds^ = f{r)df -J^- r^dnl (13) 

where dVi^ is the line element of a unit sphere. 

The effective potential for the scalar field (11) in the metric (13) is 

V.,M = /W(2^ + ^^). (14, 

The massless gauge field satisfies the equation 

V«(V„^6 - Vft^a) = 0. (15) 

After the separation of the angular variables, one can find the effective potential for the two types of the polarizations 

VGP{r)=f{r)^-^^^, (16) 

where = 1, 2, 3, . . .. 

For the fermions we do not consider the equation in the standard Schrodinger wavelike form because of the boundary 
conditions: we will calculate the grey-body factors for neutrinos and antineutrinos using the approach of [24]. That 
is why we use here the equation for the radial part in the following form 

A(r)*"(r) + ^nr) + {^-^^^^^^^^^ + " ^ = 0, (17) 
where A(r) = r^f{r) and k= 1,2,3 

III. CALCULATIONS OF THE that is, to solve the problem of classical scattering around 

ENERGY-EMISSION RATE 



In order to calculate the intensity of the Hawking ra- 
diation, one should first calculate the gray-body factors, 
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black holes with pure in-going boundary conditions at the 
event horizon. The latter is reduced to the finding of the 
S-matrix, or, simply, of the reflection or the transmission 
coefficients. 



A. Reflection coefficients 

At the event horizon we impose the boundary condition 
that corresponds to the purely ingoing wave 



$(r) cx e-'""-* cx (r - ro)"'"/^'^''''^ 



(18) 



At the spatial infinity (?■ — > oo) the two linearly indepen- 
dent solutions of the wavelike equation (7) are 

<i>(r) ~ exp(-ia;r^) -I- Zo exp(iwr*), 

where Zi and Zo are integration constants which corre- 
spond to the ingoing and outgoing waves respectively. 
Introducing the new function 

P(r) = (1 - (ro/r))'"/-'^'(''°'$(r), 

and choosing the integration constant as P{ro) = 1, we 
expand the equation (7) near the event horizon and find 
P'(ro), which completely determines initial conditions for 
the numerical integration. Then, we integrate numeri- 
cally the equation (7) from the event horizon rg until 
some distant point and find a fit for the numeri- 

cal solution far from the black hole in the form 



P{r) = Z,P,{r) + ZoPoir), 



(19) 



where the asymptotical expansions for the corresponding 
functions are found by expanding (7) for large r as 



P^{r) 
Po{r) 



1 + P^'K-^ + P^r-'- 
(l-|-P«r-i+Pi^)r-2 



The fitting procedure allows us to find the coefficients 
Zi and Zo- In order to check the accuracy of the calcu- 
lated coefficients one should increase the internal preci- 
sion of the numerical integration procedure, the value of 
R and the number of terms in the series expansion for 
Pi{r) and Po{t), making sure that the values of Zi and 
Zo do not change within the desired precision. 



B. WKB approach 

Another way to check our numerical calculations is to 
compare the fraction Zo/Zi with the result provided by 
the 6-th order WKB method, which gives quite an ac- 
curate answer for large and moderate values of I [16]. 
The WKB approach was initially used for finding quasi- 
normal modes for which it usually provides quite a good 
accuracy at the 6th order [18] 
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FIG. 3: The energy-emission rates for the Gauss-Bonnet black 
holes D = Q, a = 1/5 (blue, top line), a = 1/2 (green, middle 
line), and q = 1 (red, bottom line). 



The reflection coefficient, given by the WKB formula. 



is 



where 



\ZolZ,\^ - (1 



'2i-KK\-l 



4=6 



(20) 



(21) 



i=2 



Here Vb is the maximum of the effective potential, Vq" 
is the second derivative of the effective potential in its 
maximum with respect to the tortoise coordinate, and 
Ki are higher order WKB corrections which depend on 
up to the 2i-th order derivatives of the effective potential 
at its maximum. 

From Fig. 1 we see that the difference between the 
results found by these two methods is a fraction of a per- 
cent. Unfortunately, for some values of the parameters 
the effective potential does not have the form of a peak 
(see Fig. 2). In these cases the WKB formula cannot be 
used in the present form. That is why, despite in most 
cases the WKB formula provides a good approximation, 
we use a numerical integration procedure which works for 
any form of the effective potential. 



C. Graviton emission 

When the coefficients Zi and Zo are obtained, one can 
find the absorption probability 



|^,|2 = 1-|V^,|2 

and, then, the energy-emission rate 



dE 



duj 



exp{uj/TH) - 1 27r 



(22) 



(23) 
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FIG. 4: The energy-emission rate (left panel) and the contribution of the graviton emission of a tensor-type (right panel) for 
the Gauss-Bonnet black hole D = Q, a = 2. Blue (top) and red (bottom) lines correspond respectively to the contribution of 
1 = 2 and Z = 3 multipoles of a scalar-type (left panel) and of a tensor-type (right panel) . 
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FIG. 5: The energy-emission rates (black line) for the Gauss-Bonnet black holes D = Q, a = 1/5, 1/2, 1 (from left to right) 
together with contributions of different types of gravitons (colored lines); I = 2 scalar type (blue, top), / = 2 vector type 
(green), I = 3 scalar type (red), and I = 2 tensor type (magenta, bottom). 




FIG. 6: The energy-emission rates for D = 5 Schwarzschild (left panel) and Gauss-Bonnet black holes (a = 1/5, right 
panel). The black line is the total energy-emission rate. Scalar-type, vector-type, and tensor- type gravitons' contributions are 
red, green, and blue. The contributions of the corresponding multipole numbers are orange, yellow, and cyan. The largest 
gravitons' contribution is of a scalar-type, the smallest is of a tensor-type. 
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FIG. 7: The energy-emission rates for D = 8 Schwarzschild (left panel) and Gauss-Bonnet black holes (a = 1/5, right 
panel). The black line is the total energy-emission rate. Scalar-type, vector-type, and tensor-type gravitons' contributions are 
red, green, and blue. The contributions of the corresponding multipole numbers are orange, yellow, and cyan. The largest 
gravitons' contribution is of a scalar-type, the smallest is of a tensor-type. 




FIG. 8: The energy-emission rates for D — 9 Schwarzschild (left panel) and Gauss-Bonnet black holes (q — 1/5, right panel). 
Scalar-type, vector-type, and tensor-type gravitons' emission rates are red, green, and blue. For the Schwarzschild black hole 
the largest gravitons' contribution is of tensor-type and the smallest is of scalar-type, while for a = 1/5 Gauss-Bonnet black 
hole the largest gravitons' contribution is of a scalar-type and the smallest is of a tensor-type. The black (bottom) line is the 
energy-emission rate due to the scalar field. 



where the multiplicity factors are [15] 

(T) _ iD-l) i£+D-2)i£-l)i2£+D-3){£+D~5)l 

i V /) — 



2{£+iy.{D-3){D-5y. 

(y) _ {£ + D -3)£{2e + D -3){e + D -5)1 

(£+iy.{D-4)\ 
(S) _ i2£ + D-3){£ + D-Ay. 



N 

£\{D-3y 

for the gravitational perturbations of tensor, vector and 
scalar types. The Hawking temperature is 

/'(ro) {D - 3){24 + c^iD - ^){D - 5}} 



Th = 



47r 



Snroirl + a{D - 3){D - i)) 



D. Emission of the scalar field and the Standard 
Model fields 

The scalar field living in the bulk has the same number 
of the degrees of freedom as the gravitational perturba- 
tions of scalar type 



(SB) ^ i2£ + D-3){i + D-4y 
£l{D-3y 



For the scalar field localized on the brane the multiplicity 
factor is given 

_ (2£+lK! _ 
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FIG. 9: The energy-emission rates for D = 11 Schwarzschild (left panel) and Gauss-Bonnet black holes {a = 1/5, right 
panel). The black line is the total energy-emission rate. Scalar-type, vector-type, and tensor-type gravitons' contributions 
are red, green, and blue. The contributions of the corresponding multipole numbers are orange, yellow, and cyan. For the 
Schwarzschild black hole the largest gravitons' contribution is of a tensor-type and the smallest is of a scalar-type, while for 
a = 1/5 Gauss-Bonnet black hole the largest gravitons' contribution is of a scalar-type and the smallest is of a tensor-type. 
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FIG. 10: The energy-emission rates for D = 10 Schwarzschild (left panel) and Gauss-Bonnet black holes {a — 1/5, right panel) 
due to brane- localized fields: gauge (red, top), Dirac (green, middle), and scalar (blue, bottom). 



For the gauge field the multiplicity factor is the same 
for each polarization 



(GP) 



2e + i. 



The energy-emission rate for the neutrinos and antineu- 
trinos is the same and given by 



For the neutrinos and antineutrinos the multiplicity 
factor is 



dE 



K = l 



exp(a;/rff) + 1 2ti 



In order to calculate the absorption probability for the 
fermions we use the approach of [24] 



Note that we take into account contributions of neu- 
trinos and antineutrinos and both polarizations of the 
gauge bosons. That is why our result is two times larger 
than in [131. 
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TABLE I: Energy-emission rate for Schwarzschild (first line) and Schwarzschild-Gauss-Bonnet (a/^Q = 1/5) (second line) black 
holes. 



D 


gravitational 


scalar field 


projected SM fields 




total 


scalar 


vector 


tensor 


bulk 


projected 


s = 1/2 


s = 1 


5 


0.00133 


0.00085 


0.00041 


0.00006 


0.00105 


0.00266 


0.00464 


0.00365 




0.000322 


0.000267 


0.000055 


0.000001 


0.000230 


0.000772 


0.001259 


0.000733 


6 


0.01380 


0.00796 


0.00437 


0.00147 


0.00267 


0.01073 


0.01955 


0.01944 




0.000618 


0.000535 


0.000082 


0.000002 


0.000100 


0.001141 


0.002051 


0.001340 


7 


0.06942 


0.03229 


0.02366 


0.01347 


0.00648 


0.02972 


0.05310 


0.05935 




0.000994 


0.000833 


0.000152 


0.000008 


0.000057 


0.001832 


0.003440 


0.002528 


8 


0.27143 


0.09152 


0.09648 


0.08343 


0.01614 


0.06613 


0.11527 


0.13742 




0.001737 


0.001366 


0.000332 


0.000038 


0.000047 


0.003197 


0.006157 


0.005043 


9 


0.9975 


0.2294 


0.3440 


0.4241 


0.04216 


0.12769 


0.21768 


0.26977 




0.003490 


0.002430 


0.000897 


0.000162 


0.000050 


0.005631 


0.010781 


0.009809 


10 


3.6665 


0.5592 


1.1639 


1.9434 


0.11595 


0.22321 


0.37362 


0.47460 




0.007727 


0.004550 


0.002425 


0.000752 


0.000064 


0.009633 


0.018190 


0.017976 


11 


13.7745 


1.4696 


3.8936 


8.4113 


0.34326 


0.36257 


0.59793 


0.77210 




0.018552 


0.008939 


0.006678 


0.002935 


0.000094 


0.015795 


0.029331 


0.030885 



IV. RESULTS 

In Fig. 3 one can see that the energy-emission rate 
per unit frequency w for gravitons strongly decreases as 
a/r^ increases, leading to suppression of about five times 
for a/rg = 1. Examples of contributions of various mul- 
tipoles and of various types of gravitons (scalar, vector 
and tensor) are given in Figs. 4, 5, 6, 7, 8, 9. These 
are given for £) = 6 as an example, while from the table 
I, we can see the relative contributions for all types of 
particles for various numbers of space-time dimensions. 
In Fig. 10 one can see contributions of brane-localized 
Standard model fields into the energy-emission rate: for 
D > 6 the larger spin of a field, the larger is the corre- 
sponding energy-emission rate per unit frequency around 
its maximum. 

An important observation for our future discussion is 
about the role of the tensor type of gravitons in the 
evaporation process. For vanishing GB coupling the ten- 
sor type of gravitons corresponds to the highest energy- 
emission rate, which leads to the dominance of tensorial 
gravitons among all particles in the Hawking radiation at 
large D (see Figs. 8, 9). For non zero a this is not true 
anymore: the energy-emission rates of tensorial gravitons 
become strongly suppressed and the Standard Model par- 
ticles dominate in the evaporation at high D. 

Finally, let us justify here in more detail why we chose 
the units of the fixed event horizon and not the units of 
the fixed black hole mass. 

First, we would like to mention that from the values 
in units of the black hole radius, using a simple formula 
which we shall show in the next section, we are always 
able to recalculate values of the energy-emission rate in 
units of the black hole mass or Planck mass. The inverse 
formula is not so simple because the mass appears in 
the left-hand side as well. Second, in order to calculate 
the energy-emission rate in units of the mass, we must 
fix in some way the real Planck mass, which we do not 
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FIG. 11: The energy-emission rate as a function of mass M 
for D = 10 Schwarzschild (top line) and a = O.lTeV"^ GB 
(bottom line) black holes, M« — ITeV. According to the 
recent experiments M, > 0.6TeV (for D = 10) [14]. For 
larger M, the plots are qualitatively the same. 

yet know. Third, the purpose of the table I is to show 
that, independently on the new Planck mass, the energy- 
emission rate decreases some orders. This is a general 
result, which we would lose if we fix the Planck mass in 
some way. In order to show this effect we present a figure 
in the units of TeVs to illustrate how the emission rate 
depends on the back hole mass (Fig. 11). 



V. DISCUSSIONS 

From Fig. 12 and Table I one can see that the energy- 
emission rate decreases quickly (in fact, almost exponen- 
tially) with a . In order to understand why this happens 
for relatively small corrections to the black hole geome- 
try let us use the WKB approach and consider the cases 
of small and large a; of separately. 

The numerical data, which we obtained here (Table I) , 
do not use the small uj expansion of [13] and are, thereby. 
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FIG. 12: Energy-emission rate of the tensor-type gravitons as 
a function of a for D = 6. 
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FIG. 15: Grey-body factor of the tensor-type gravitational 
perturbations as a function of q {D = 6, £ — 2, ujro — 1.5). 
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FIG. 13: Expansion of the tensor-type perturbation potential 
D = 5,£ = 2. 



more accurate. 

When oj^ < Vmax, the energy-emission rate decreases 
due to the suppression of the grey-body factor for ten- 
sor gravitons. In order to show this we use the WKB 
approximation. For <C Vmax the WKB formula reads 



exp 



(^-2 J dr,^V{r,)-uj^^ , (24) 



where integration is performed between the two turning 



points V{r^,) 



On Fig. 13 we see that the leading 



order of a increases the height of the potential barrier, as 
well as the distance between the turning points. Then, 




0.05 



0.10 



0.15 



0.20 



FIG. 14: Black hole temperature as a function of a for D = 
11. 



from (24) it is clear that that the grey-body factor de- 
creases exponentially as the effective potential for tensor 
gravitons grows. 

In the eikonal approximation for cj^ < Vmax one has 



1 + g-2'r(V„„-a.2)/^_2y4'„ 



As the height of the potential increases with growing 
a, the dominant contribution comes from the numera- 
tor and the grey-body factor for tensor gravitons is 

\A\^ oc e"^(")", C{lj) > 0. 

For large values of w it is evident that \A\'^ « 1. Thus, 
the dominant contribution to the energy-emission rate at 
large to comes from the temperature, which decreases as 
(see Fig. 14) 



Th = Th{0)(1 



(D - l)(7:>-4)a 



2r2 



Therefore the energy-emission rate decreases exponen- 
tially with a 



dtduj 



Ecu 
2n 



dtdoj 



xe 

a=0 



u{D-l)iD-4:)a/{2r^TH) 



Let us compare the Hawking radiation of Schwarzschild 
and GB black holes produced due to particle collisions of 
the same energy. The relation between the black hole 
mass M and its radius ro is given by 



{D - 3)iD - A)a 



2r2 



(roAf,)^-3(i?-2)^(^-3)/2' 

(25) 

where M, is the true fundamental scale of gravity. In 
units of the black hole horizon the energy-emission rate 
^ ^q'^ . If we measure the black hole radius in units of 
TeV^^, the energy-emission rate is measured in units of 
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TeV"^. In order to convert the energy-emission rate in 
the units of TeV/s, we divide it by the Plank constant 
h = 6.58 • IQ-'^^TeV ■ s. For a = O.IAC^ ^nd D = 10, 
the total energy-emission rate of GB black holes can be 
as much as 10^ times smaller (see Fig. 11) than that of 
the Schwarzschild one. For other values of D the increase 
in the lifetime of black holes can be easily calculated in 
the same way. 

Thus, we have shown that even small GB corrections 
to the £)-dimensional Schwarzschild geometry lead to the 
great increasing of the lifetime of black holes, up to quite 
a few orders. This is certainly not enough for accreting of 
matter and thus is not dangerous for experiments at the 
LHC, yet, it will produce a potentially observable effect. 

A natural question is, what is the fraction of the contri- 
bution of the grey-body factor into suppression of emis- 
sion of the tensorial gravitons, and consequently into 
overall lifetime of black holes. An exact answer is, how- 
ever, not as easy as it drastically depends on a number 
of conditions: the mass of the black hole M, the value 
of the GB coupling a, the number of space-time dimen- 
sions D and is strongly connected with the dominance of 
a particular type of particles for fixed a, M , D. 

A brief hint could be given by the plot of the grey- 
body factor as a function of a, which is shown in Fig. 
15. For example, the grey-body factor for D = 6 can be 
suppressed at the rate of about 35% for sufficiently small 
a/rQ 0.1, which is one tenth of its threshold values 
of instability [20] . Going beyond very small values of a 
would require inclusion of corrections of higher orders of 
curvature in the action. Thus, the first factor of the two 
which was mentioned in the abstract of this article, i.e. 
the quick cooling of a black hole when a grows, is defi- 
nitely a dominant factor of suppression of the Hawking 
radiation. 



Let us note that the suppression of the energy-emission 
rate at high w due to the cooling of a black hole was ob- 
served in [11] for the Standard Model particles. Though 
the largest suppression of the graviton emission reported 
in [11] is at M < M,. In this regime however one cannot 
trust the solution (2), which is gravitationally unstable 
[20]. The other factor, the decreasing of the grey-body 
factors of the tensorial gravitons, seemingly increased the 
suppression of the graviton emission at small a. 



VI. CONCLUSIONS 

We have shown that the widely accepted approxima- 
tion of higher dimensional black holes by their classical 
Schwarzschild- Tangherlini model is not good, when one 
considers the Hawking radiation around a black hole. In- 
tensive Hawking emission of gravitons, as well as of other 
particles, is suppressed by many orders, when one takes 
into consideration small quantum Gauss-Bonnet correc- 
tions. Consequently, the lifetime of quantum corrected 
black holes is many orders larger than it is expected ac- 
cording to the current literature [25] . This makes further 
investigations of Hawking radiation of higher curvature 
corrected black holes appealing. 



Acknowledgments 

R. K. was supported by the Alexander von Hum- 
boldt Foundation, Germany. A. Z. was supported by 
Fundagdo de Amparo d Pesquisa do Estado de Sao Paulo 
(FAPESP), Brazil. 



[1] P. Horava and E. Witten Nucl. Phys. B 460 506 (1996). 

[2] I. Antoniadis, Phys. Lett. B 246, 377 (1990); N. Arkani- 
Hamed, D. Dimopoulos and G. Dvali G Phys. Lett. B 
429 263 (1998); Phys. Rev. D 59 086004 (1999); I. Anto- 
niadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, 
Phys. Lett. B 436 (1998) 257 [arXiv:hep-ph/9804398]; 
L. Randall and R. Sundrum Phys. Rev. Lett. 83 3370 
(1999); Phys. Rev. Lett. 83 4690 (1999). 

[3] D. G. Boulware and S. Deser, Phys. Rev. Lett. 55, 2656 
(1985). 

[4] R. Konoplya, Phys. Rev. D 71, 024038 (2005) 
[arXiv:hep-th/0410057]; E. Abdalla, R. A. Konoplya and 
C. Molina, Phys. Rev. D 72, 084006 (2005) [arXiv:hep- 
th/0507100]; A. Zhidcnko, Phys. Rev. D 78, 024007 
(2008) [arXiv:0802.2262 [gr-qc]]. 

[5] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[6] E. Jung and D. K. Park, Nucl. Phys. B 731, 171 (2005). 
Mod. Phys. Lett. A 22, 1635 (2007) [hep-th/0612043]. 

[7] S. Chen, B. Wang, R. K. Su and W. Y. Hwang, JHEP 
0803, 019 (2008) [arXiv:0711.3599 [hep-th]]. 

[8] J. A. Prost, J. R. Gaunt, M. O. P. Sampaio, M. Casals, 



S. R. Dolan, M. A. Parker and B. R. Webber, 

arXiv:0904.0979 [hep-ph]. 
[9] D. Ida, K. y. Oda and S. C. Park, Phys. Rev. D 71, 

124039 (2005); Phys. Rev. D 73, 124022 (2006). 
[10] H. Nomura, S. Yoshida, M. Tanabe and K. i. Maeda, 

Prog. Theor. Phys. 114, 707 (2005). 
[11] T. G. Rizzo, Class. Quant. Grav. 23, 4263 (2006) 

[arXiv;hcp-ph/0601029]. 
[12] J. Grain, A. Barrau and P. Kanti, Phys. Rev. D 72, 

104016 (2005) [arXiv:hep-th/0509128]; 
[13] V. Cardoso, M. Cavaglia and L. Gualtieri, Phys. Rev. 

Lett. 96, 071301 (2006) [Erratum-ibid. 96, 219902 

(2006)] [arXiv:hcp-th/0512002]. 
[14] P. Kanti, Lcct. Notes Phys. 769, 387 (2009) 

[arXiv:0802.2218 [hep-th]]. 
[15] P. Kanti, H. Kodama, R. A. Konoplya, N. Pappas 

and A. Zhidenko, Phys. Rev. D 80, 084016 (2009) 

[arXiv:0906.3845 [hep-th]]. 
[16] B. F. Schutz and C. M. Will Astrophys. J. Lett 291 L33 

(1985); S. Iyer and C. M. Will Phys. Rev. D 35 3621 

(1987); R. A. Konoplya, Phys. Rev. D 68, 024018 (2003) 



13 



[arXiv:gr-qc/0303052]; J. Phys. Stud. 8, 93 (2004). 

[17] R. A. Konoplya and A. Zhidcnko, Phys. Rev. D 81, 
124036 (2010) [arXiv: 1004. 1284 [hcp-th]]; Phys. Lett. B 
686, 199 (2010) [arXiv:0909.2138 [hep-th]]. 

[18] O. P. F. Piedra and J. de Oliveira, arXiv: 1006.3802 [gr- 
qc]; R. A. Konoplya, Phys. Lett. B 679, 499 (2009) 
[arXiv:0905.1523 [hcp-th]]'; M. 1. Liu, H. y. Liu and 
Y. X. Gui, Class. Quant. Grav. 25, 105001 (2008) 
[arXiv:0806.2716 [gr-qc]]; J. F. Chang, J. Huang and 
Y. G. Shen, Int. J. Theor. Phys. 46, 2617 (2007); 
Y. Zhang and Y. X. Gui, Class. Quant. Grav. 23, 6141 
(2006) [arXiv:gr-qc/0612009]; H. T. Cho, A. S. Cor- 
nell, J. Doukas and W. Naylor, arXiv:0912.2740 [gr-qc]; 
R. A. Konoplya, K. Murata, J. Soda and A. Zhidenko, 
Phys. Rev. D 78, 084012 (2008) [arXiv:0807.1897 [hep- 
th]]. 

[19] A. Ishibashi and H. Kodania, Prog. Theor. Phys. 
110, 901 (2003) [arXiv:hep-th/0305185]; H. Ko- 
dama and A. Ishibashi, Prog. Theor. Phys. Ill, 29 

(2004) [arXiv:hcp-th/0308128]; [arXiv:gr-qc/0312012]; 
R. A. Konoplya and A. Zhidcnko, Phys. Rev. D 78, 
104017 (2008) [arXiv:0809.2048 [hcp-th]]; Phys. Rev. 
Lett. 103, 161101 (2009) [arXiv:0809.2822 [hep-th]]; 
Nucl. Phys. B 777, 182 (2007) [arXiv:hep-th/0703231]. 

[20] G. Dotti and R. J. Glciscr, Phys. Rev. D 72, 044018 

(2005) [arXiv:gr-qc/0503117]; Phys. Rev. D 72, 124002 



(2005) [arXiv:gr-qc/0510069]; 

[21] T. Takahashi and J. Soda, [arXiv:1008.1618 [gr-qc]]; 
Phys. Rev. D 80, 104021 (2009) [arXiv:0907.0556 [gr- 
qc]]. 

[22] R. A. Konoplya and A. Zhidenko, Phys. Rev. D 77, 

104004 (2008) [arXiv:0802.0267 [hep-th]]. 
[23] P. Kanti and J. March-Russell, Phys. Rev. D 66, 024023 

(2002) [arXiv:hcp-ph/0203223]; Phys. Rev. D 67, 104019 

(2003) [arXiv:hcp-ph/0212199]; P. Kanti, Int. J. Mod. 
Phys. A 19, 4899 (2004) [arXiv:hep-ph/0402168]. 

[24] M. Casals, S. R. Dolan, P. Kanti and E. Winstanley, 
JHEP 0703, 019 (2007) [arXiv:hep-th/0608193]. 

[25] F. Moura and R. Schiappa, Class. Quant. Grav. 24, 361 
(2007) [arXiv:hep-th/0605001]; S. Creek, O. Efthimiou, 
P. Kanti and K. Tamvakis, Phys. Rev. D 76, 104013 
(2007) [arXiv:0707.1768 [hep-th]]; Phys. Lett. B 635, 39 

(2006) [arXiv:hep-th/0601126]; V. Cardoso, M. Cavaglia 
and L. Gualtieri, JHEP 0602, 021 (2006); [arXivdicp- 
th/0512116]; T. G. Rizzo, JHEP 0506, 079 (2005) 
[arXiv:hep-ph/0503163]; J. L. Hewett, B. Lillic and 
T. G. Rizzo, Phys. Rev. Lett. 95, 261603 (2005) 
[arXiv:hcp-ph/0503178]; A. Barrau, J. Grain and 
S. O. Alexeyev, Phys. Lett. B 584, 114 (2004) [arXivihep- 
ph/0311238]. 



